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Non-dissipative Josephson current through nanoscale superconducting constrictions is carried by 
spectroscopically sharp energy states, so-called Andreev bound states. Although theoretically pre- 
dicted almost 40 years ago, no direct spectroscopic evidence of these Andreev bound states exists 
to date. We propose a novel type of spectroscopy based on embedding a superconducting constric- 
tion, formed by a single-level molecule junction, in a microwave QED cavity environment. In the 
electron-dressed cavity spectrum we find a polariton excitation at twice the Andreev bound state 
energy, and a superconducting-phase dependent ac Stark shift of the cavity frequency. Dispersive 
measurement of this frequency shift can be used for Andreev bound state spectroscopy. 



Supercurrents through mesoscopic or nanosized 
Josephson junctions are mainly carried by spectroscop- 
ically sharp subgap bound states, so-called Andreev 
bound states (ABS)[P, as predicted theoretically in 
Refs. [2J [3J. When the weak link is a quantum dot with 
a single or a few discrete levels supercurrents are also 
predicted to flow mainly through ABS [I]. ABS come 
in pairs, one state above and one below the Fermi level. 
The two ABS of the pair have opposite dispersion with 
the superconducting phase difference over the junction 
and carry supercurrent in opposite directions across the 
junction. These two states form a well-defined two-level 
system that has been suggested as a qubit [5] [5]. Be- 
fore such a qubit can be realized, experimental detection 
and characterization of this engineered two-level system 
should be carried out. To our best knowledge, how- 
ever, no one has to date reported experiments with direct 
spectroscopic proof for the existence of Josephson-current 
carrying ABS [7]. 

A single- wall carbon nanotube (swCNT) embedded be- 
tween two superconducting metal leads can support a 
supercurrent [3 M EH HH E3 E3J. The electronic en- 
ergy levels of the swCNT, formed through size quanti- 
zation, can be tuned in and out of resonance with the 
lead Fermi levels by gating the swCNT. Such supercon- 
ducting swCNT transistors [TTJ can further be switched 
from a Coulomb blockade regime, to a Kondo regime, to 
a weakly interacting Fabry-Perot regime by changing a 
back gate voltage [12]. The potential for applications 
of such swCNT quantum-level junctions was demon- 
strated through the fabrication and detailed functional 
control of a nano-SQUID, involving two gated swCNT 
junctions with controlled on-off states as well as con- 
trolled 0-7r SQUID states [13]. Similar control has been 
demonstrated using semiconducting nanowire junctions 

[HCElirB]. 

We propose a method for direct ABS spectroscopy, 
based on dispersive measurement of a polaritonic state 
formed by the ABS strongly coupled to a QED cavity 
mode. Consider two superconducting leads connected by 



a molecule with one resonant level, as shown in Figure 1. 
This setup can be realized in the swCNT transistor by 
tuning the voltage of the gate electrode [TT] [12] [13] . The 
superconducting proximity effect leads to a split of the 
resonant level into a pair of ABS, where the ABS level 
splitting can be controlled by tuning the superconduct- 
ing phase difference in a SQUID setup. When the gate 
electrode is coupled to an LC-oscillator, induced quan- 
tum fluctuations of the gate potential leads to dressed 
Andreev bound states. By tuning the ABS into reso- 
nance with the cavity, and measuring the shift of the 
cavity base-frequency, the ABS energy can be observed. 
The complete ABS energy dispersion with respect to the 
superconducting phase can be extracted by using a cav- 
ity with a variable base-frequency [T7 1 ITS ] f!9 | [20 ] [21] . 
The suggested dispersive method does not involve real in- 
terlevel transitions and could be simpler to realize than 
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FIG. 1: Superconducting density of states for tp = 0.67T re- 
solved into right, N+(e), and left, N-(e), current-carrying 
branches. The lightly shaded background is the density of 
states of the resonant level when the contacts are in the nor- 
mal state. The electron-oscillator coupling leads to dressed 
Andreev bound states with sidebands separated in energy by 
the oscillator frequency uj . The inset shows a sketch of the 
considered system. 
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methods based on microwave absorption [22 . Another 
advantage of the suggested spectroscopy is that it makes 
it possible to directly access and characterize individual 
ABS without measuring dc contact current- voltage char- 
acteristics [23 . Similar dispersive measurement method 
has earlier been applied to detect quantum states of su- 
perconducting charge qubits |24j . 

The Hamiltonian for the coupled electron-oscillator 
system depicted in the inset in Figure [l] is 

H = Hl + Hr + Hl eve l + Ht + H osc + Hi eve i- osc (1) 

and will be quantified in terms of creation operators 
for the three excitations of the system: reservoir elec- 
trons, ct (reservoir a = L, R, momentum k, and 

spin a); dot electrons, dy, and the oscillator mode w. 
The two superconducting reservoirs are described by a 
standard BCS Hamiltonian H a = Y,ka^c\ a a c ka , a + 

~^2 k ^A Q c^,| cfi^-ki a + h.c^j , where h.c. denotes hermi- 
tian conjugate, is the quasiparticle dispersion, and 
A Q = A(T)e' lVa are the order parameters of the reser- 
voirs that have the same temperature dependent gap 
A(T), but have a tunable superconducting phase differ- 
ence if = ipn — ipi between them. The reservoirs are 
coupled by tunneling through a single non-interacting 
molecular level described by Hi eV ei — J2a £ dd] J d G , where 
the level energy 6d(V g ) is tunable by the gate volt- 
age. The lead-to-level tunneling is described by, H T — 
J2ka,a { V k^acl a a d a + h.c)j . We consider symmetric, 
spin-independent coupling between the level and the two 

leads, Vka-.L = v^a.R = v5k,k F , which corresponds to 

i i 2 

an effective tunneling rate at each barrier T = -^-Afp, 
where Nf is the normal state density of states of the leads 
at the Fermi-level. The oscillator mode with frequency lo q 
is described by the Hamiltonian H osc — uu wb. The cou- 
pling of the molecular level to the oscillator is described 
by a linear interaction iJ iet , ei _ osc = X(b + tf)dld a , 
where A is the coupling strength. We will assume 
throughout this paper that the molecular level is aligned 
with the Fermi levels of the leads, i.e. Ed = 0. 

We solve the Hamiltonian (TXT) treating the tunneling to 
infinite order in the hopping [25 , 26, 27 and the electron- 
oscillator coupling perturbatively in a self-consistent 
Born-approximation [2SJ CSS 113 ED GUIMI ■ The problem 
reduces to the coupled molecular level/oscillator system. 
Superconductivity modifies the electronic state on the 
molecular- level by splitting it into two branches (s = ±1) 
each described by a retarded Green's function 



G R {s) 



Q(e) 



z R (e)e R + sA cos(y>/2) ' 



(2) 



where the energy-renormalisation factor z s {e) is defined 



i R 



by z R (e)e 

v/|A| 2 - (e R ) 2 /2T. The state of the oscillator is given 



ft(e) (e R - £*,(e)) , with Q(e) 



by the retarded Green's function 

2<Jo 

(w R ) 2 -LJ 2 -2io U R (cu)' 



D R {u) 



(3) 



The electron-oscillator coupling enters via the retarded 
self-energies and Ii R (uf) which are determined 

self-consistently by numerical iteration. T, R (e) and 
H R (cu) are functionals of Gf(e) and D R (u) [and of 
Gf(e) = -2iImGf(e)n e (e), with n e (e) = tanh(e/2T) 
and D k (uj) — —2ilmD R (oj)nb(oj), with nb{uS) — 
coth(w/2T), Gf = (Gf )*] and defined as 

Ef( £ ) =z^[D K (Lu)oGf(e - W )+D»oGf( £ - W )] , (4) 
H^)^zA 2 ^[G fl s (e)oGf( £ -c) + G^( £ )oG^ £ -4) 

s=±l 



where the convolution is defined as a(x) o b(x — y) = 
^r a ( x )H x ~ v)- Up° n reaching self-consistency the 
Josephson current is calculated as 



S = ±l 



f —s Im 




Loo 27T 


[ n(s) \ 



n e (e). (6) 



Inelastic coupling to the environment is introduced 



phcnomcnologically above by e 



~R 



i"f and u> 



to + in. The parameter 7 <C A describes a residual 
phase-breaking scattering rate in the superconducting 
reservoirs. In absence of an electron-oscillator coupling it 
is 7 that limits the life time of the ABS. K describes the 
finite (long) life time of the mode originating from a finite 
(but high) quality factor, Q, such that k — luq/Q -C u)q. 
In the present calculations we set 7 = n = 10 _3 A. 

Without coupling to the oscillator, the Josephson ef- 
fect in this system is well known 0]. The Green's 
function-amplitudes, G R , describe two different quasipar- 
ticle branches that form on the dot, one having an ABS 
below (s = +1) and one having an ABS above (s = —1) 
the Fermi level, see Figure 1. The two branches con- 
tribute to the Josephson current ([6| in opposite manner. 
For sub-gap energies Gf (\e\ < A) contributes to the cur- 
rent in the positive direction for s — 1 (left to right over 
the junction) and vice versa for the s = — 1 branch. The 
continuum part of G R (\e\ > A) contributes to the current 
in the opposite direction compared with its correspond- 
ing sub-gap part. 

The electron-oscillator interaction couples the two 
quasiparticle branches. This is seen explicitly in the 
energy-renormalization factor z s of one branch which is 
modified by the self-energy S^ s of the other branch. In 
Figure [2] we show a fully self-consistent calculation of the 
Andreev spectrum as function of superconducting phase- 
difference. For the chosen parameters there is a resonance 
between the oscillator and the ABS, i.e. 2Ea{<p) = UJ 0) 
at cp w 0.97T. Away from resonance, the ABS is shifted 



E A ^E A = E A - 



\2E A 



2E A 



(7) 
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FIG. 2: (a) The spectrum of Andreev-bound states as 
function of superconducting phase-difference [tp) for electron- 
oscillator coupling A = 0.1A and oscillator frequency u) = 
0.2A. The tunneling rate is F = A and the temperature is 
T — 0.05A. The side bands to each ABS are due to the 
dressing of the ABS by the electron-oscillator coupling. The 
sidebands of the ABS with branch index s = ±1 belong to 
the quaasiparticle branch with index s =p 1. (b) The current- 
phase relation for different electron-oscillator couplings for 
T — 0.05A. (c) The reduction of current due to coupling to 
the oscillator as a function of phase for the same temperature 
and electron-oscillator couplings as in panel b). 

as compared with the case without coupling to the oscil- 
lator, but with retained phase-dependent spectral weight 
Ae of the state at e = ±Ea(<p)- The effective electron- 
oscillator coupling in the sub-gap region is given by the 
product XAe- In equation ([T]) thermal occupation factors 
enter in the combinations $± = nf,(w ) ±n e (£U). Apart 
from the shifted ABS we find satellite resonances at e = 
sEa^lOo with spectral weig hts \X 2 A e $± s /{2sE a ±uj ) 2 . 
It is only at resonance, IE a — oj , the satellite on ei- 
ther branch with spectral weight cx interfere with the 
main ABS of the same branch index. This is seen as a 
precursor of an avoided crossing in Figure [2] 

The Josephson current-phase relation is presented in 
panels b and c of Figure [2] for the case that the elec- 
tron and oscillator systems are in thermal equilibrium. 
Due to the different magnitudes of the thermal factors 
<I> ± ($+ « 2 and w for T < E A ,w ), the current 
contribution of the satellites is dominated by the satellite 
with weight cx This satellite reduces the Josephson 
current by ~ X 2 A E /(2E A + oj ) 2 , which is of the order of 
a few percent of the full current for our parameter values. 
There is no dramatic signature in the current-phase re- 
lation of an emerging anticrossing at resonance because 
interference occurs between states carrying current in the 
same direction, and moreover, these states have spectral 
weights shared between them drawn from the original 
ABS, and the population of the states is largely phase- 
independent. 

The Josephson current through the level modifies the 
oscillator spectrum giving it a phase-dependence that is 
shown in Figure |3ja) for the same parameters as for the 



FIG. 3: (a) The spectra of the oscillator plotted as function 
of superconducting phase-difference with the same parame- 
ters as in Figure [2] A mode emerges at uja(<p) = 2E A (f) and 
develops to a polariton when in resonance with the oscillator 
(inset, 0.8tt < <p < n in steps of 0.0125tt). The dashed box 
indicate where the spectra in the inset are taken. Away from 
resonance there is a Stark shift 5u}(ip) of the base frequency, 
(b) The Stark shift is shown at different electron-oscillator 
coupling strengths for T = 0.05A and (c) at different temper- 
atures for a electron-oscillator coupling strength A = 0.05A. 
In both cases uj — 0.2A. The dashed line in panels b)-c) 
indicate the resonant phase. 

ABS spectrum in Figure [2] The polarization H R (oj) gives 
the possible collective excitations that are supported by 
the electronic system. Iterating the self-consistency equa- 
tions once we get 



n R (u) = 2A 2 



A\ 



4E/ 



-,n e {E/ 



(8) 



for the retarded phonon self-energy. As may be expected 
there is a mode with the phase-dispersion loa = 2Ea(^>) 
originating from the transitions between the two ABS 
and the subsequent emission or absorption of the energy 
w A - A spectral weight ~2\ 2 A E [<ko' 2 /(4E\-u 2 ) 2 }n e (E A ) 
for this collective excitation is vanishingly small away 
from resonance. Near resonance the oscillator mode and 
the excitation interact strongly and an avoided crossing 
appears with a frequency split at 2Ea = lo q 




2X 2 A 2 E n e (E A ) 



(9) 



where S = 2Ea — u is the detuning. The analytic 
estimate of the split, 2\Ae\[2mJJ^a) ~ 0.096, for a 
nearly phase-independent Ae ~ 0.34 when r = A, is 
in good agreement with the self-consistently determined 
split, w 0.1 A, extracted from the inset in Figure]^ a) at 
f = 0.9 7T. It is important to note that this split is signifi- 
cantly larger than the intrinsic broadening of the spectral 
features (~ k) and signals a strong-coupling regime in the 
sense of cavity-QED, i.e. XAe ^> k. 

Away from resonance, the electron-oscillator coupling 
gives a negative phase-dependent Stark shift of the oscil- 
lator base frequency. We plot this shift in Figure [3Fb)-(c) 
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both as a function of electron-oscillator coupling strength 
A for T = 0.05A and as a function of temperature for 
A = 0.05A. Our analytic estimate in equation ^ gives 
in the limit of large detuning 



iv (ip) = u + 5u((p) = uj — 2 



X 2 A 2 E n e (E A ) 
IE A - u 



(10) 



This analytical expression for 5ui(tp) is a good approxima- 
tion when the continuum contributions can be neglected, 
i.e. for r/A <C 1, while in general one must use nu- 
merics to extract duj(ip). Dispersive measurement of the 
acquired phase-dependent resonance frequency of the os- 
cillator, uj (ip), gives a possibility to detect the position 
of the polariton resonance and hence define the energy of 
the ABS. Furthermore, by sweeping the base frequency 
of the oscillator the phase dispersion of the ABS energy 
can be detected. 

Assume that we bridge two aluminium superconduc- 
tors by a gated swCNT, as shown in the inset of Figure [l] 
The gate is part of an LC-circuit which leads to an oscil- 
lation of the gate voltage with frequency uj = 1/ \f~LC. A 
practical value for the resonator frequency is ~ 10 GHz, 
for which the polariton resonance is located well inside 
the superconducting gap (for aluminium Aai ~ 50 GHz), 
and at the same time the ABS energy splitting is large 
compared to the temperature below lOOmK. In our cal- 
culations we neglect electron-electron interactions. This 
can be done if T is large compared to A, which can ex- 
perimentally be realized by tuning T by a back gate [12 
to approach the weakly interacting Fabry-Perot regime. 
The main modification of the one-iteration approxima- 
tion in the case T 3> A is that Ae becomes strongly phase 
dependent, tending to (A/2T) sin(y>/2) as T/A grows. 
This gives the effective coupling XAe ~ AA/r close to 
the resonance. This coupling must be large compared to 
the intrinsic oscillator damping, AA/r » k = lo /Q, 
in order to resolve Sui(tp). The strength of the bare 
coupling A is determined by the capacitive interaction 
between the gate and the dot, and it is proportional 
to the ratio of corresponding capacitances, C g /Cs, and 
can be expressed through the oscillator frequency as 
A = (C g /C s )(S c /8S L ) 1 / 4 o; (E c and E L are charging 
and inductive energies of the oscillator, respectively). For 
swCNT superconducting contacts the gate capacitance 
can be comparable with the capacitances of the contacts 
to the leads, having values of tens of aF [52]. Thus the 
coupling A can be on the order of 10% of the oscillator 
frequency as assumed in our calculation. When the gate 
of the contact is connected to a superconducting cavity, 
the quality factor of the oscillator can be of order 1000 
- 10000 [HI El El 120] • Given our calculated value for 
effective ABS-oscillator coupling, XAe, we get the ratio 
XAe/k ~ 100 — 1000, indicating that the strong coupling 
regime is indeed feasible, the resolution of the proposed 
spectroscopy should be very favorable to encourage ex- 
periments. 
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